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On Jacobi’s Method of Facilitating the Numerical Solution of 
Equations arising in the Theory of Secular Perturbations. 

By H. C. Plummer, M.A. 


1. The equations to be considered here are of the type 

A = o 


where A is a symmetric determinant which may be written in the 
form 


An A, A I2 , . . , A ln 
A 2 i, A 2 2 A, ... A 2n 


A 




A 


n 2) 



and which is no other than the familiar discriminating determi¬ 
nant which occurs in so many branches of mathematics. As it 
arises in the determination of the mean motions of the perihelion 
and the node, the number of rows and columns corresponds to 
the number of the planets. In the application to the solar 
system the determinant is of the eighth order ; and little re¬ 
flexion is necessary to see that its direct reduction would entail 
enormous labour. Hence Jacobi’s method * of approximating 
to the numerical values of the roots without reducing the deter¬ 
minant is of great practical interest. 

2. The principle of Jacobi’s method may be briefly recalled. 
The vanishing of A is the condition of the consistency of the 
system of n linear equations 

(An \)x t * 4 “ A I2 x 2 -f- • • • d~ A In x n = o 

A„i x 1 -f- A n2 x 2 -f . . . -f (A )in — h)x tl = o 

where' A ra = A sr . Let [i] and \Jc\ denote the ith and kth equations. 
The substitution 

»< = Xi cos a ~ Xk sin a 
x k = Xi sin a + Xk cos a 

is made and [il and [&] are replaced by [i] cos a + \Jc\ sin a and 
— [i] sin a + [&] cos a. Then if 

cot 2a = (A a - A u )/2A a 

the eliminant of the equations thus transformed is a symmetric 
determinant of exactly the same form as before. If A' rs be the 
constituent corresponding to A rs in the transformed determinant, 
we have 


* * fiber ein leichtes Verfahren die in der Theorie der Sacnlarstorungen 
vorkommenden Gleichungen numerisch aufzulosen/ WerJce , vol. vii. p. in ; 
Crelle, yo], xxx. (1846), p. 64. 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjoumals.org/ at University of Ulster at Coleraine on April 30, 2015 






1902MNRAS..63...90P 


91 


Dec. 1902. the Numerical Solution of Equations etc . 


A-'ik — A! Hi — 0 

A'ii 2 + A'** 2 = A,? + A kk 2 + 2A ik 2 
AV + AV = A ? f + A^ 2 , jjiijik. 

The result of the transformation is therefore to annihilate two 
equal constituents and to leave unaltered the sum of the squares 
of the other constituents outside the leading diagonal. The 
process can be repeated until the sum of the squares of the con¬ 
stituents becomes very small except in the leading diagonal, the 
constituents in which approximate to the factors of the determi¬ 
nant, and will accordingly give the roots of the equation. The 
approximations to the roots after a certain number of transforma¬ 
tions can be completed by a different process, also due to Jacobi. 

3. It may be interesting to study the theory underlying this 
method in a more general manner. For this purpose let the 
determinant A be multiplied by 



ct 12 

... 

» 

Cl 2 i 

a 2 2 

. . . 

« 2n 


G'm n 2 • • • a n n 

the mode of composition being by rows of each determinant in 
such a way that rows of A appear in the corresponding rows of 
8A. The constituent in thejth row and sth column of the product 

1 = n 

— ^ A-ji a 8 i 

i — I 

Next let the determinant 8A be multiplied by 8, the mode of 
composition being by columns of 8A and rows of 8. The con¬ 
stituent in the r th row and s th column of the final product 

j = n r-i *= n —j 

== ^ a r j 2 A ji cigi X. a s j 
j = i Li = i J 

i = n j = n j—n 

= 2 2 A a H a rj —X. 2 a rj a v -. 

i = I j = I j = X 

This expression is unaltered by the interchange of r and s, and 
therefore the resulting determinant is symmetrical. 

4. The determinant 8 2 A is not only symmetrical but is also of 
precisely the same form as A, provided that 

Yv = i,(f = i, 2.») 

3 = 1 

2 a r j a sj = o, (V, s = 1, 2, r zfi s) 

3 = 1 

That is, provided that h is the modulus of an orthogonal trans¬ 
formation scheme 

x x — cCh y~i 4" a 2 x y% H - • • • -h a nJ y n 


x n — cf>j n y 1 4 " a 2n y% -H • • • a nn y n 
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which is such that 


!2 +^ 2 2 +... + xj = yi 2 + yS +... + yj 

If these conditions are fulfilled, with their well-known conse¬ 
quences, and if 

i—nj—n 

2 2 A^ a 3i a rj = B fS 

i — 1 j — l 


then S 2 A = 


Bn—A, B I2 , ... B ln 
B 2I , B 22 —X, ... B 3ft 



B„ 


B„ 


— \ 


is identically equal to A, for S 2 = i. 

5. A number of identities may be deduced from this result. 
The coefficient of (—A) 7i_I in A is 

Ah+A 22 + • • • +A nn =:$A rr 

Hence 2 A rj . = SB n . 


Again the coefficient of (—\) n 3 

= S(A r ,A ss ~A, a 2 ) 

Hence 

[3 A,.,.] 2 — 2 ^(A rr A ss '—A r , 2 ) 

= [SB rr ] 2 —2 S(B n .B a *—B r8 2 ) 

which gives 

2 A )T 2 -m2 A,/ = ^B r , 2 + 2^B,. S 2 

That is to say, the sum of the squares of the constituents is the 
same in both determinants. Moreover, the sum of the squares of 
the constituents in the rth row or column of 8 2 A is 


2B rs 2 = 2 (S 2 A n a si a rj ) 2 

s 3 j i 

= SSSSSA^A kl a rj a rk a si a s l 
s j i Jc l 

i j k 

for %a ti a a = o unless l = i, in which case the sum is unity. 

s 

6. The conditions which the n 2 coefficients of an orthogonal 
transformation must satisfy are -}n(n + 1) in number. They can 
be so chosen as to satisfy in addition the \n(n— 1) conditions. 

i—n j—n 

B,.„ = 2 2 Ajj a si a rj = o, (r, s = 1,2, . . .,n; r zfi s) 

i = 1 j= 1 


Then all the constituents outside the leading diagonal of 
vanish, and the roots of the equation in \ become B n , B 22 , . 


B, w , where 


i=n j~n 

B r) . = S % Aj l a ri a r 


1 — 1 ? — 1 


S 2 A 
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Bat the complete solution of the equation A = o is required in 
order to find the coefficients of the substitution. The additional 
relations given above show that 

A II a; I 2 -fA 22 a? 2 2 + . . . + 2A I2 x I x 2 + . . . EE . . . + B nn y n 2 

i.e. the substitution is one which transforms the quadratic form 
into a sum of square terms. Hence 

(A ii —X)^ i 2 + (A 23 — \)x 2 2 + . . . + 2A l2 x 1 x 2 

+ • * * ■ X)y J 2 + . . . + ^)Vn 

If now X = B rr , y r disappears, and therefore 

,,[(A„-X)| + a „| + . . . +A„| ; ] 

■ ■ • +x -%] 

+ . . . — o 

Hence 

(An -|- A I2 a r2 + • • . "f" A ln a rn == o 

N 21 a n + (A 22 —X(a r2 + . . . A 2n a rjl = o 


These lead again to the equation A = o. The successive substitu¬ 
tion of the roots of the latter equation in X in the set of linear 
equations will yield the n sets of n coefficients. 

7. It is thus seen that the complete transformation requires a 
knowledge of the roots of the very equation whose solution is to 
be effected. But it is possible to use a partial transformation 
instead of a complete one. The set of quantities x 19 ... ,x n can be 
divided into two sets x 19 ... 9 x m and x m+I9 . . . , x n of which the 
former is transformed as before, while the latter remains 

unchanged. This means that the terms of %%Aj^Xj which 

i j 

involve the first set only are transformed into a sum of squares, 
those involving the second set only are unchanged, and the 
product terms involving one quantity from each set retain their 
form with changed coefficients. The transformation as a whole, 
as well as in part, remains orthogonal, but the coefficients of the 
substitution are such that a si = o (s zfi i; s or i>m) and 
cOgg = i (s > m). Three cases must be considered separately, and 
these results are easily deduced from the formula of § 3 : 

(1) If r>?n 9 s>m , then 

Bj -,3 = A rs 

(2) If r> m, s^j>m, then 

i — m 

B r8 = % Ayi a si 

i = 1 

(3) If s^>m 9 then 

j — m i = m 

i>rs =: ^ ^ A-jfl'sfl'rj 

j- 1 i- 1 
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8. It is now possible to choose the m 2 quantities a si (s3>m, i^>m) 
according to the method of § 6 so as to make the terms zfi s) 
of § 7 (3) vanish. In this way all the constituents of the minor 
formed of the first m rows and columns of A outside the leading 
diagonal can be made to vanish. Since by interchange of rows 
and columns the constituents of the leading diagonal can be 
made to take any order, it is clear that any minor on the leading 
diagonal can be transformed in this way. If the minor is a 
determinant of the mth order, the solution of an equation of the 
mth degree is required for the transformation and m(m— 1) con¬ 
stituents can be made to vanish. To the transformed minor the 
results of §§ 3-6 apply. The complementary minor on the 
leading diagonal will remain unaltered. As regards the out¬ 
standing constituents it is easily deduced either from § 7 (2) or 
from the formula in § 5 that the sum of the squares of con¬ 
stituents lying in a line parallel to a row or column of the 
transformed minor remains unaltered. This is the generalisation 
of the fact on which the principle of Jacobi’s method depends. 
The sum of the constant terms in the axis of the transformed 
minor is unchanged, while the sum of their squares is increased 
by the sum of the vanishing constituents. Outside this minor 
the sum of the squares of the constituents is left, line for line, 
unchanged. Hence repetition of the process must cause the axial 
terms to approximate to the factors of the determinant. 

9. Jacobi’s method amounts to the use of the case in which 
m = 2. The auxiliary angle a which is introduced in § 2 has of 
course a simple geometrical meaning. It is the angle between 
the coordinate axes and the axes of the conic 

Ai,x° + 2 A ik xy +A kk y* = i 

and its connexion with the substitution theory is obvious. In 
Jacobi’s method it is always possible to remove the numerically 
greatest pair of constituents outside the diagonal at each trans¬ 
formation. Hence if S is the sum of the n(n— 1) squares after the 
first transformation the square of the greatest >S j(n —2)(w+i), 
and each further transformation will reduce the sum of the 
squares faster than in the ratio 

2 

1 ; 1—- v --- 

(n—2)(n-\- 1) 

In the general case of m substitutions the m(m —1) constituent 
to be removed cannot be chosen independently, and after a certain 
number of transformations it may be necessary to include terms 
previously reduced. Short of this event the ratio of decrease of 
the sum of the squares may be expected to be at least 

i j m(m — 1) 

(n—?n)(n + m —1) 
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Doubtless Jacobi's case is the simplest, but it would be 
theoretically possible to make m = 4, a case involving the 
solution of a biquadratic and removing twelve constituents. The 
case of m = 3, depending on a cubic and removing 6 constituents, 
would seem to be quite practicable and useful at least in the 
earlier transformations. In any case, when the factors of any 
axial minor are known, all the constituents of that minor outside 
the axis can be removed. The necessary substitutions might 
probably be facilitated by the use of an arithmometer. But 
after the work of Harzer, who employed Jacobi’s method, it is 
unlikely that an occasion for employing such processes will 
recur. * 

10. The principle of Jacobi’s method contains implicitly a 
proof that the roots of a discriminating determinant are real. 
The proof supposes the method repeated an infinite number of 
times. The idea of a limit may be avoided by the foregoing 
theory. For this it is only necessary to notice that if the extra¬ 
diagonal terms of a minor of any order be removed as before 
explained, it is possible to extend the process to a minor of an 
order higher by unity. That is to say, the process can be 
extended successively from a minor of the second order, for 
which it is always possible, to the complete determinant. For 
let the extra-diagonal constituents of A be zero except in the 
first row and column. Then 

_ -A. I2 2 -^-I3 2 

(A„-X) (A 22 -X) (A ix —\)(A 33 —X) • • • 

where II is the product of the diagonal terms. In A put X = A n .. 
Then 

A = A Ir 2 (A 22 A rr ) (A 33 A rr ) . . . (N nn A rr ) 

Suppose that A 22 , A 33 . . . are in increasing order of magni¬ 
tude. If now X is put equal to —00, A 22 , A 33 , . . ., A nn , +00 
in turn, it is easily seen that the resulting values of A are of 
alternating sign, and therefore A has real roots situated between 
the values, assigned to X. And if A 22 , A 33 , . . . are not all 
different, the roots are still real. For if, for example, A 22 = A 33 , 
the expanded form of A shows that A 22 is a root, and the order 
of A can be depressed by writing A 22 —X only once in the 
diagonal and \/[A I2 2 -f A I3 2 ] in the corresponding places of the 
first row and column. Hence, generally, a determinant of 
the form considered has real roots, and its extra-diagonal con¬ 
stituents can all be removed by a real substitution. Thus by 
successive steps it is proved that the roots of any discriminating 
determinant are real. It is not difficult to see that a similar 
method could be applied to the more general case of Langrange’s 
determinant. 

University Observatory, Oxford: 

1902 December 6. 
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Ephemeris for Physical Observations of 
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